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ABSTRACT 

We present Ks, H & J-band photometry of the very highly irradiated hot Jupiter WASP-12b using 
the Wide-field Infrared Camera on the Canada-France-Hawaii telescope. Our photometry brackets 
the secondary eclipse of WASP-12b in the Ks and H-bands, and in J-band starts in mid-eclipse 
and continues until well after the end of the eclipse. We detect its thermal emission in all three 
near-infrared bands. Our secondary eclipse depths are 0.309to;Qi2% m Ks-band (24c), 0.176to;Q2i% 
in H-band (9a) and 0.13lto'o29% m J-band (4<r). All three secondary eclipses are best-fit with a 
consistent phase, (f>, that is compatible with a circular orbit: 0=O.4998^q'qoo7- The limits on the 
eccentricity, e, and argument of periastron, uj, of this planet from our photometry alone are thus 
|ecosw| < 0.0040. By combining our secondary eclipse times with others published in the literature, 
as well as the radial velocity and transit timing data for this system, we show that there is no 
evidence that WASP-12b is precessing at a detectable rate, and show that its orbital eccentricity is 
likely zero. Our thermal emission measurements also allow us to constrain the characteristics of the 
planet's atmosphere; our Ks-band eclipse depth argues strongly in favour of inefficient day to nightside 
redistribution of heat and a low Bond albedo for this very highly irradiated hot Jupiter. The J and 
H-band brightness temperatures are slightly cooler than the Ks-band brightness temperature, and 
thus hint at the possibility of a modest temperature inversion deep in the atmosphere of WASP-12b; 
the high pressure, deep atmospheric layers probed by our J and H-band observations are likely more 
homogenized than the higher altitude layer probed by our Ks-band observations. Lastly, our best- 
fit Ks-band eclipse has a marginally longer duration than would otherwise be expected; this may be 
tentative evidence for material being tidally stripped from the planet - as was predicted for this system 
by Li & collaborators, and for which observational confirmation was recently arguably provided by 
Fossati & collaborators. 

Subject headings: planetary systems . stars: individual: WASP-12 . techniques: photometric- eclipses 
- infrared: planetary systems 



1. INTRODUCTION 

Multiwavelength constraints on the thermal emission 
of hot Jupiters are crucial to precisely defining the spec- 
tral energy distributions of these planets and under- 
standing their energy budgets. Interestingly most hot 
Jupiter thermal emission detections to date have not 
been at the blackbody peaks of these planets, but at 
longer wavele ngths with the Spitzer Space Telescope (A 
> 3 u rn; e.g. iCharbonneau et al~l 120051 iDeming et al. I 
2005). Probing shorter near-infrared wavelengths at 
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the blackbody peaks of these planets has only re- 
cently been proven feasible first through space-based 
observations with t he Hubble Space Telescope (HST; 
i Swain et al. Il2009af), and then from the ground (e.g . 
de Mooii fc Snellen I l200l iSing fc Lopez-Morales 1 120001 : 
ICillon et al. II2009D . Our program to detect near-infrared 
thermal emission from the hottest of the hot Jupiters has 
also been successful using the Wide-field Infrared Cam- 
era (WIRCam) on the Canada-France-Hawaii Telescope 
(CFHT) to detect the Ks-b and thermal emission of: 
TrES-2b dCroll et al Il2010ai). TrES-3 b including an H- 
band upper- limit ( Croll et al. I l2010bh . and two eclipses 
of WASP-3b. includin g a limit on its temporal variability 
(jCroll et al. Ilin prep.f ). 

In the near-infrared multiple band detections have only 
been performed on a handful of occasions; such multiple- 
band detections were performed in narrow wavelength 
regimes from space via spectroscopy with H ST for HD 
209458 and HD 189733 ([Swain et al. 2009a ; b), and ar- 
guably recently from the gro und for HD 189733 using the 
Infrared Telescope Facility (|Swain et al. Il20l0l ). as well 
as from the ground using the Very Large Telescope in 
the H & K- bands for the highly irradiated hot Jupite r 
WASP-19b (|Anderson et al. ll2010HGibson et al. Il2010t ). 
Multiple band detections in the near-infrared are there- 
fore rare compared to the frequent multiple -band detec- 
tions at longer wavelengths using the IRAC (iFazio et al. I 
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2004), IRS (|Houck et al. Il200l . or MIPS (IRieke et al. I 
20041) instruments on the Spitzer Space Telescope. 



Multiwavelength thermal emission measurements with 
Spitzer have revealed a wealth of information, includ- 
ing that the most highly irradiated exoplanets seem 
to ha rbour hot stratospheres and temperature inver- 
sions (|Knutson et al. I 2008at iCharbonneau et al. 1120051 : 
iMachalek et al. I 120081 Knutson et al. I l2008bl ). One 
could imagine that obtaining multiwavelength con- 
straints on a planet's thermal emission in the near- 
infrared could be equally informative. Furthermore the 
near-infrared is also an ideal place to directly constrain 
these planets' pressure-temperature profiles at depth, 
dayside bolometric luminosities and the fraction of the 
incident stellar radiation that is transported from the 
tidally locke d day to nights ides deep in these planets' 
atmospheres (Barman 2008). 

Here we continue our program using WIRCam on 
CFHT to detect thermal emission from some of the 
hottest of the hot Jupiters. Our target was the highly 
irradiated hot Jupiter WASP-12b. The discovery of the 
inflated, transiting exoplanet WASP-12b was of immedi- 
ate interest to those attempting to measure the loss in 
flux during the secondary eclipses of hot Jupiters in the 
near-infrared - this was because WASP- 12b circles a late 
F-typ e star with a period of only ~26 hours ()Hebb et al. I 
2009). It is thus exposed to extremely high stellar in- 
solation, with an incident flux of ^9xl0 9 ergs~ 1 cm~ 2 . 
The planet is also one of the most inflated hot Jupiters, 
with a radius of Rp^l.&Rj and a favourable plane t-to- 
star radius ratio (i? P /i?»-0.12; IHebb et al. "1 120091 ) . It 
should be heated to an equilibrium temperature of over 
^2500 K assuming isotropic reradiation and a zero Bond 
albedo^. For these reasons it was predicted to display 
near-infrared thermal emission on the order of 0.1-0.3% 
of the stellar flux in the J, H & Ks near-infrared bands, 
assuming isotropic reradiati on and a zero Bond albedo. 
ILopez-Morales et al. I (|2010D have already reported a de- 
tection of the secondary eclip se of WASP-12b in z' -band 
(0.9 urn), and more recently iCampo et al. I (|2010l ) have 
presented detections of two eclipses in the f our IRAC 
channels for WASP-12b. ICampo et al. I (|2010), however, 
did not report the eclipse depth s for WASP- 12b, and 
for re asons discussed below the ILopez-Morales et al. I 
(2010) detection has recently been called into question. 
Thus the atmospheric characteristics of WASP- 12b re- 
main largely unconstrained. 

In addition to receiving extremely high stellar insola- 
tion, WASP-12b is intriguing because the combination 
of its close proximity to its s tar and its putativ e original 
eccentricity (e=0.049±0.015: lHebb et al. 1120091 ) suggests 
that it could be precessing at a rate that is detectable 
with current instruments. Such a putative pr e cessio n 
signal was recently claimed by ICampo et al. I f2 010). 
Altho ugh the IRAC eclipses reported by ICampo et al. 1 
(2010) suggest an ecosw constraint similar to that 
expected for a circul a r orbi t (ecosw=-0.0054±0.0030), 
ILopez-Morales et al. I (|2010f ) had earlier reported an 
eclipse detection that was considerably offset from a cir- 
cular orbit (|ecosa; |=0.016t£^g). While at first glance 
the two measurements are inconsistent, if the planet pre- 

8 The Bond albedo is the fraction of the bolometric flux reflected 
from the planet compared to the incident bolometric radiation. 



cesses this is not the c ase. By combining thei r seco ndary 
eclipses with those of ILopez-Morales et al. I (|2010l ). to- 
gethe r with the original radial velocity data for the sys- 
tem (jHebb et al. I [2001 . as well as a series of transit- 
time measurements from the original detection paper 
and groun d-based amateurs (from the Exoplanet Transit 
Database; iPoddanv et~aT~ll2010D . ICampo et al~l (|2010D 
show that a precessing orbital model best-fits the data 
with 2(T confidence. The authors caution that this de- 
tection is heavily dependent on the sec o ndary eclipse 
offset reported by ILopez-Morales et al. I (|2010f ) . Even 
more recently, radial velocity observations of WASP- 
12b have suggested that the eccentricity of WASP-12b 
is small (e=0.017j:o.gio ; IHusnoo ef^IH 120101) and likely 
zero, constraining the ICampo et al. I (I2010D precession 
signal and calling into question the ILopez-Morales et al. I 
(2010) eclipse detection. Nevertheless, the best-fit ec- 
centricity of WASP-12b remains non-zero, and thus this 
planet could be prec essing at a much slower rate than 
ICampo et al. I (|2010D claim. The definitive nail in the 
coffin on the claim that WASP- 12b is precessing at a 
detectable rate, will thus only result from further detec- 
tions of this planet's secondary eclipse well seperated in 
time from the original eclipse detections. 

Also, recently preliminary evidence was presented that 
material from WASP-12b may be being tidally stripped 
from the planet an d may pos s ibly f orm a circumstellar 
disk in this system. iLi et al.l 1)20101 ) predicted this sys- 
tem may have such a disk from material overfilling the 
Roche lobe of WASP-12b, because WASP-12b's observ ed 
radius in the optical (R p ~l.79 Rr, IHebb et al. II2009D is 
al ready close to its 2.36 Rj Roche lobe radius (as quoted 
in IFossati et al. Il2010at ). That WASP-12b may exhibit 
material overfilling it Roche lobe and a circumstellar disk 
from this material has recently received possible confir- 
mation from HST Cosmic Origins Spectrograph (COS) 
observations of this s ystem. From these observations 
iFossati et al. I (|2010aD find increased transit depths in 
the ultraviolet when compared to the optical, indicative 
of material surrounding WASP- 12b overfilling its Roche 
lobe and blocking out a larger fraction of the stellar flux 
at these wavelengths. In addition they observe an early 
ingress of th e transit of WASP- 12b in their near ultra- 
violet data; IFossati et al. I (|2010aD intrepret this early 
ingress as a putative sign of previously stripped mate- 
rial from WASP-12b forming a circumstellar disk. These 
putative signs of a disk are interesting to ob servers in 
the ne ar-infrared, specifically the K-band, as ILi et al.l 
(2010) predicted that such a disk in this system may 
exhibit CO emission as bright as 10 mJy at 2.292 ixm. 
WASP- 12 does not, however, display a significant near- 
infrared excess (jFossati et al. Il2~010b[ ). 

Here we present detections of WASP-12b's thermal 
emission in the Ks (24er), H (9<r) and J-bands (4c). 
Our J-band detection is the first thermal emission mea- 
surement in this band. Our photometry favours a cir- 
cular orbit for WASP-12b (ecosw=-0.0007^;^^). By 
combining our secon d ary e clipse times with those o f 
ILopez-Morales et al. I (|2010l ) and lCampo et al. I feOlO). 
as well as the radia l veloc ity data of IHebb et aL I pool) 
and IHusnoo et al~l (|2010l) . and all the transit-time data 
for the system, we are able to show that not only is there 
no evidence to date that WASP- 12b is precessing at a 
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Figure 1. The normalized flux from our target star and refer- 
ence stars for our Ks-band photometry (top two panels), our ri- 
band photometry (middle two panels), and our J-band photometry 
(bottom panels). For each set of panels the top panel displays the 
flux from the target star (black) and the reference stars (various 
colours) that are used to calibrate the flux of WASP-12b in the var- 
ious sets of photometry. The bottom panels in each set of panels 
displays the residuals from the normalized flux of the target star 
corrected by the normalized flux of the reference stars. 

detectable rate, but also that the orbit of WASP- 12b is 
likely circular. Our analysis also allows us to constrain 
the characteristics of the atmospere of WASP-12b; our 
Ks-band eclipse depth argues in favour of inefficient re- 
distribution of heat from the day to nightside, while our 
J and H-band observations seem to be probing deeper, 
higher pressure atmospheric layers that are slightly more 
homogenized. We also show that our Ks-band photome- 
try may feature a longer than expected eclipse duration 
that could arguably be interpreted as evidence for mate- 
rial streaming from the planet or a circumstellar disk in 
this system. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained observations with WIRCam on CFHT of 
WASP-12 (J-10.48, iJ-10.23, if-10.19) on 2009 De- 
cember 26, 27 and 28 in the J, H & Ks-bands respectively. 
Our J-band observations on Dec. 26 lasted for 3.9 hours 
and started in mid-eclipse and persisted for 2.2 hours af- 
ter the end of eclipse. Our observations on Dec. 27 and 
28 lasted for 6.0 hours in H-band and 6.2 hours in Ks- 
band, respectively, evenly bracketing the predicted sec- 
ondary eclipse of WASP- 12b. Numerous reference stars 
were also observed in the 21x21 arcmin field-of-view of 
WIRCam. The telescope was defocused for our various 
observations to approximately 1.5mm (J-band), 1.8mm 
(H-band), and 2.0mm (Ks-band), resulting in the flux 
of our target star being spread over a ring ~19, ^23, 
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Figure 2. The root-mean-square of our out-of-eclipse photome- 
try (solid line) following the corrections documented in J2]f or °ur 
Ks-band photometry (top), our H-band photometry (middle), and 
our J-band photometry (bottom). The dashed line in each panel 
displays the one over the square-root of the bin-size expectation 
for gaussian noise. 

and ^26 pixels in diameter (6, 7 and 8") on our ar- 
ray. For each observation, as the telescope temperature 
changed over the course of the night, we used the fo- 
cus stage model and kept the defocus amount constant, 
thus achieved a stable PSF over the entire observation 
set. We used "Staring Mode" for our J and Ks-band 
observations where we do not dither for the duration of 
our observations; for the H-band eclipse the queue obser- 
vations mistakenly used micro-dithering which featured 
small 0.5 pixel shifts between consecutive exposures. The 
exposure times for our J, H & Ks-band observations were 
5-seconds. The effective duty cycle after accounting for 
readout and for saving exposures was 34%. 

For both observations the data was reduced and aper- 
ture photometry was performed o n our target s tar and 
our reference stars as discuss ed in iCroll et al. I (12010H) 
(with the details provided in ICroll et al. Il2010al) . We 
used an aperture with a radius of 17 pixels for our Ks- 
band photometry, and 16.5 pixels for our H and J-band 
photometry. We used an annulus to define the sky with 
an inner radius of 22, and an outer radius of 34 pixels 
for all our photometry. We ensured that these choices of 
aperture were optimal by testing smaller and larger aper- 
ture sizes in increments of 0.5 pixels and ensuring these 
choices displayed the smallest root-mean-square (RMS) 
outside of occultation and the least time-correlated red- 
noise. 

Following our aperture photometry we correct the flux 
of our target st ar with a number of n earby reference stars 
as discussed in ICroll et al. I (|2010af) . We used 22, 7, and 
17 reference stars to correct our J, H and Ks-band eclipse 
photometry, respectively. The normalized flux of WASP- 
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12 and the various reference stars that are used to correct 
the flux of our target star are displayed in Figure [TJ For 
our Ks-band photometry we corrected our photometry 
for a small trend with the x, and y pixel position of the 
target star on the chip 0. We didn't notice such trends 
in our H and J-band photometry. 

For our H-band photometry the airmass, X, was high 
at the start of the observations (X~1.9), and fell to 
X^1.2 by mid-eclipse. We noticed a downward trend 
in our H-band photometry following the correction with 
nearby reference stars that appeared to be correlated 
with airmass. We found this effect was reduced, but not 
removed, for our H-band photometry by correcting the 
flux of WASP- 12 with reference stars solely on the same 
WIRCam chip as WASP-12; this downward trend in flux 
of our target star compared to the reference stars is still 
apparent at the start of our H-band photometry. To re- 
duce the impact of these systematic data we scale-up the 
errors of the first ^25 minutes of data for our H-band 
photometry by a factor of 1.3. 

The root-mean-square (RMS) of our photometry per 
minute following the above corrections improved from 
lO.OxKT 3 to 0.86xl(T 3 in Ks-band, 3.8xl(T 3 to 
0.90X10" 3 in H-band and 1.93xl(T 3 to 0.69xl(T 3 in 
J-band. To evaluate the impact of systematics and the 
presence of red-noise in our photometry we bin our data 
and compare the out-of-eclipse photometric precision to 
the gaussian noise expectation of one over the square- 
root of the bin-size (Figure [2]). The Ks-band data bins 
down very close to the gaussian noise limit, while the 
J-band data bins down marginally above this limit; the 
H-band data is worse, possibly due to the systematics 
introduced by the micro-dithering. To ensure we do not 
underestimate the uncertainties in our model parameters 
we employ the lWinn et al. I (|2008f ) method to account for 
time-correlated red-noise in our photometry. We scale- 
up the uncertainties on our individual data-points by a 
factor fj; (3 is equal to the factor that the binned out-of- 
eclipse RMS scales above the gaussian noise expectation 
in the absence of red-noise. We use a binning time of 
~12 minutes. For our H-band photometry we exclude 
the hrst 15 minutes of data from this calculation, due to 
the obvious systematic that we believe to be correlated 
with airmass that does not appear to affect the rest of the 
photometry For our photometry /3 is equal to 1.7, 1.3 
and 1.1 for our J, H & Ks-band data, respectively. We 
note that our observations are still well above the pre- 
dicted photon noise RMS limit of 2.7xl(T 4 , 2.2 xl(T 4 
and 3.1 xlCP 4 per minute in the J, H & Ks-bands, re- 
spectively. 

3. ANALYSIS 

For many of our other WIRCam data-sets we have ob- 
served residual background trends in the reduced data 
that seems to a ffect our target star s differently than our 
reference stars (jCroll et al. l l2010aHb1 . lin prep I ) . For our 
Ks, H & J-band photometry these backgrounds, B/, dis- 
played a near-linear slope. We fit our Ks, H & J-band 
data-sets with linear backgrounds of the form: 



Phase 

0.4 0.42 0.44 0.46 0.48 0.5 0.52 0.54 0.56 0.58 0.6 0.62 0.64 



B f = 1 + ci +c 2 dt 

The correction is described in lCroll et al. I Ij2010al ) 
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Figure 3. CFHT/WIRCam photometry of the secondary eclipse 
of WASP-12b observed in the Ks-band on 28 December 2009. The 
top panel shows the unbinned lightcurve with the best-fit secondary 
eclipse and background from our MCMC analysis of the Ks-band 
data with the fixed eclipse duration (red line). The second panel 
shows the lightcurve with the data binned every ~7.0 minutes and 
again our best-fit eclipse and background. The third panel shows 
the binned data after the subtraction of the best- fit background, 
Bf, along with the best-fit eclipse model. The bottom panel shows 
the binned residuals from the best-fit model. 
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Figure 4. The same as figure [3] except that the data is our ri- 
band photometry obtained on 27 December 2009. 

where dt is the interval from the beginning of the ob- 
servations and c\, C2 and C3 are fit parameters. We use 
Markov Chain Monte Carlo (MCMC) fitting to fit for 
our background as well as a second a ry ecl ipse model 
calculated from the IMandel fc Agol I (|2002t ) algorithm 
without limb darkening. We fit for the background, the 
depth of the secondary eclipse (AF) and the offset that 
the eclipse occurs later than the expected eclipse cen- 
ter (t ff set)- Our Markov Cha in Monte Ca r lo meth od is 
discussed in lCroITI (f2006h and lCroll et aTl (I2010afl . We 
obtain our stellar and planetary parameters for WASP- 
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Table 1 

Best-fit secondary eclipse parameters 



Parameter Ks-band H-band J-band Joint Ks-band MCMC Joint MCMC 

MCMC MCMC MCMC MCMC variable eclipse variable eclipse 

Solution Solution Solution Solution duration solution duration solution 



reduced \ 2 
AF Ks 
AF H 
AFj 

toff set (min) a 

teclipseKs (B JD-2450000) 
teMpseH (B JD-2450000) 
tecUpseJ (BJD-2450000) 
ClKs 
C2Ks (d' 1 ) 

C\H 
C2H (d- 1 ) 

ClJ 

C2J (d- 1 ) 



11/1 

<&H (hours) 
Tbks (K) 
Tbh (K) 
Tbj (K) 
ecos(oj) a 
e sin(oj) 

fKs 
f H 
fj 



n 7 o 9 +0.003 
U.IOZ_ 001 

o.3iol°;°il% 

n/a 
n/a 
-1.3__ 
5194.9351 
n/a 
n/a 
0.00116 
0.001 
n/a 
n/a 
n/a 
n/a 
0.4994 
n/a 
2.93 
2988 
n/a 
n/a 
-0.0013 



0.435 



+ 1.5 
1.2 

+0.0010 
-0.0008 



0.00014 
0.00014 
+0.001 
0.001 



+0.0009 
0.0008 



j+45 
-46 



0.482 



+0.0015 
0.0015 
n/a 
+0.030 
0.029 
n/a 

n/a 



+0.002 
0.001 
n/a 

0-l80±g;S|96 



0.358 



2.0 1 



n/a 

1+2.1 
-2.6 

5193.84*61+°;°^ 
n/a 
n/a 

n/a 

(-0. 00022 



0.00021 
■0.001 



0.00109] 

-0.003; _ 

n/a 

n/a 

u.ouio_ 0017 

n/a 

2.93 

n/a 

2765j 

n/a 

0020+ 0021 
u.uuzu_ 0021 

n/a 

n/a 

ot-o+0.037 
u.ooo_ 036 

n/a 



0.003 
0.002 
n/a 

0.126+H?°% 



-2.4^1; 

n/a 
n/a 

5192.7515 
n/a 
n/a 
n/a 
n/a 

0.00157 
-0.014 



+0.0031 



+0.00027 

-0.00035 

+0.003 

0.003 

n 4Q87+ - 0028 
0.498 r_ 0028 

n/a 

2.93 

n/a 

n/a 

28331 

0024+ 0044 
-u.uuz4_ 0044 

n/a 

n/a 

n/a 

0.389 



ft ttqq^u.uuz 

u.ooo_ Q01 
0.3091°°^% 

0.1761°;°^% 



0.131 



+0.027 



-0.029 ' 

-0 71 1 ' 3 

5194.9356+° 1 ;°°°? 

5193.84411°;°°°° 

5192.75271°;°°°° 

0.00116+°;°°°" 

001+ ' 001 
u.uui_ 001 

00104+"""" 2 ' 

U - UUiU4 -0. 00024 

o no3+°-° 01 

00163+""""''" 

-0.0l5l°;°°i 
0.49981°;°°°? 

n/a 

2.93 



+0.091 
0.087 



2985l 4 ° 
2748l 74 
28601^ 
-0.00071°;°°" 

0.480l°;°p 

0-3451°;°^ 
0.405l°-°f? 



n 797 +0.003 
u - ,zl -0.001 

0.3Hl°-°»% 



n/a 
n/a 



-0 9+ 1,4 
5194.9355+°;°° 1 1 ° 
n/a 
n/a 

00126+°'°°° 16 

U - UUiZD -0. 00014 

001 + 001 
u.uui_ 001 

n/a 

n/a 

n/a 

n/a 

0.4997; 

1.109 
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a Wc account for the increased light travel-time in the system ILoeb 120051) , and use the best-fit period for the non-precessing case reported by lCampo et al. I 
f2010h . 
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Figure 5. The same as figure [3] except that the data is our J- 
band photometry obtained on 26 December 2009. Note that the 
photometry is a partical eclipse only, and starts in eclipse and 
extends well out of eclipse. 



12 from lHebb et al. I (|2009). whil e the planetary period 
and ephemeris are obtained from iCampo et al. I (|2010D 
from their non-precessing best-fit. 



The best-fit secondary eclipes from our individual 
MCMC analyses with a fixed eclipse duration are pre- 
sented in Figures [3l [4] and [5] and the best-fit eclipse 
parameters are presented in Table [TJ along with associ- 
ated parameters, such as the best-fit phase, <j>, and the 
barycentric julian date of the eclipse center in the ter- 
restrial time formatQ, t ec n pse . The phase dependence 
of these fits are presented in Figure [5] We also perform 
a joint analysis of the three secondary eclipses with a 
common offset from the eclipse center (t // se t); the fit 
parameters are thus AFr-s, AFjj, AFj, i //set, and ci, 
and C2 in each band. The resulting best-fit parameters 
of this joint fit are listed in Table [TJ 

We also repeat our fit for our Ks-band photometry, 
our highest signal-to-noise photometry, and for the joint 
analysis while fitting for an additional parameter - the 
duration of the secondary eclipse, We parameterize 
this by the duration of the eclipse divided by the dura- 
tion of the transit, $////, using the duration of the tran- 



sit ($/~2.93 h) reported bv iHebb et al~l pQQjj ). The 
results from this fit are presented in Table Q] We do not 
fit our J-band or H-band data individually with this ad- 
ditional parameter, as the J-band data is a partial 
eclipse and thus the duration of the secondary eclipse is 
degenerate with an offset of the eclipse center, and the 

10 As calculated using the routines of Eastman e t al. I Q2010T ). 
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Figure 6. The 68.3% (ltx; solid-line), 95.5% (2a; dashed-line), and 99.7% (3ct; dotted-line) credible regions from our individual MCMC 
analyses with a fixed eclipse duration of our Ks-band photometry (left), H-band photometry (middle), and J-band photometry (right). The 
"x" in the middle of the plots marks the best-fit point from our MCMC analyses. 



H-band data suffers from additional time-correlated sys- 
tematics that could lead to erroneous conclusions. 

4. DISCUSSION 

We strongly detect all three secondary eclipses in the 
three near-infrared bands that we observed in. The in- 
dividual analyses of our three eclipses confirm that all 
three secondary eclipses are fit with a consistent phase 
(Table [1]); thus the best-fit parameters from our joint 
analysis are similar to the parameters returned by the 
analyses of the individual eclipses. We therefore quote 
the results of the joint analysis below. The best-fit eclipse 
depths from our joint analysis is 0.309lp' r U 9% in Ks- 
band, 0.176^;^% in H-band and 0.131 
band. 



0.012' 
+0.027iv 
-0.029 70 



4.1. Eccentricity and Precession of WASP-12b 
The best-fit phase of the joint analysis 



-0 4998+ ' 0008 
-u.ia»o_ .0007- 



centncity, e, 
ecosw=-0.0007 



is 

The resulting limit on the ec- 
and argument of periastron, u>, is 
a result that is consistent with 



+0.0013 



-0.0013' 



a circular orbit and the iCampo et al. I (|2010D re- 
sults. This value is incon sistent, however, with the 
ILopez-Morales et al. I (120101) ecoso; result . The discrep 



ancy betwe en the Lopez-Mora l es et al~l (|2010f ) result 
and that of ICampo et al. I (|2010f ) and our own could be 
due to WASP- 12b precessing - we explore this possibility 

b elow. 

ICampo et al. I (|2010l ) performed an analysis of the re- 
ported transit times and secondary eclipse times and pre- 
sented tentative evidence that WASP-12b may be pre- 
cessing at an observable rate, u> = 0.02 ± 0.01° d , 
with a period as short as 40 years. The primary evi- 
dence for the preces sion was the ground-based se condary 
eclipse detection of ILopez-Morales et al. I (|2010f ). which 
occured late by approximately ~15 minutes (at a phase 
of 0=O.51OO^;^6i using the lHebb et al~|[2009l ephemeris 
and period). 

We repeat the ICampo et al. 1 (|2010h precession analy- 
sis adding in ou r three secondary ecli pse detections. We 
summarize the ICampo et aT~l (20101) precession model 
that we employ here. The mid-transit time of the N th 
transit, Tjv, in our precessing model is predicted to occur 
at: 

e p a 

Tyv = T Q + P S N (cosujn — coslOq). (2) 

T a and uj are the transit time and argument of peri- 
astron at the reference epoch, wjv is the argument of 



25 
20 
15 
10 

5 


-5 
-10 - 



WASP-12b 



3000 



3500 



4000 4500 
BJD-2450000 



5000 



5500 



25 
20 
15 
10 
5 

-5 
-10 



WASP-12b 



3500 



4000 4500 
BJD-2450000 



5000 



5500 



Figure 7. Transit (black points) and eclipse (red points) 
times for WASP-12b compared to the best-fit orbital models 
for the precessing case (top), and the non-precessing case (bot- 
tom). The best-fit models for the transit times (dotted black 
line) and the eclipse times (solid red line) are also shown. Both 
diagrams show the observed-minus-calculated (O-C) times from 
a linear ephemeris calculated using T and P\ the secondary 
eclipse O-C times are compared to a calculated eclipse centre of 
T + v- T ne eclipse points in the two panels are from left- 
to-right, the ICampo et"aT~l H2010T) Spitzer / IRAC eclipses (BJD- 
2450000~4750). the ILopez-Morales et al. I (t20Tot) eclipse (BJD- 
2450000^5000), while the last three red-points are the eclipse times 
reported here (BJD-2450000~5200). 

periastron of the TV*' 1 transit, P s is the sidereal pe- 
riod, P a is the period between successive periastron pas- 
sages, and e has already been defined as the eccentric- 
ity. P a is not an independent variable, but is related 
to the sidereal period, P s , and the constant precession 
rate, Co: P a = y ^ . The argument of periastron 

of the N th transit is simply ujn = u>(Tn — T ) + oj . 
Equation [2] is solved iteratively for Tjv after it is ex- 
panded to fifth order in e (as shown in equation (22) of 
iRagozzine fc W olf 2009]). We fit the radial vel o city d ata 
from iHebb et al. I (|2009l) and iHusnoo et al~l (I^ IOM 
and t he transits listed in Table 2 of ICampo et al~] 
(2010) as well as four additional, recent transits n 



11 IHusnoo et"aT~l f2 010) argu e that there may be correlated red 
noise in the lHebbetal - 1 | | 2009) radial- velocity data, possibly due 
to a systematic offset in the RV zero- point from night to n ight. As 

~T l(200gh data by 



a result we scale-up the e rrors for the 



a factor of 8 and those of Husnoo et a 



Hebb et al. . 
I H20101) by a factor of 2 to 



account for possi ble offsets betwee n these two data-sets. We refer 
the reader to [Husnoo et al. (2010) for further discussion. 

12 The additional transits have mid-transit times 
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Table 2 

WASP- 12b orbital parameters 



Parameter 


Precessing 


Non- P recessing 


Non- Precessing Case without 




Case 


Case 


the Lopez-Morales et al. (2010) eclipse 



P (days) 
e 

T (BJD-2450000) 

(°) 
U) (° d" 1 ) 



X 2 
BIC 



+0.0000004 
0.0000004 
+0.01365 
0.00063 
+0.0001 
0.0002 
184.5 a 
4.5 

+0.002 
0.001 
0000+ 0003 

u.uuuu_ .o 00g 
-0.004+ 016 



1.0914239 
0.00095 
4508.9769 
-90.9 
0.003 



-0.013 

85.5 
106.6 



1.0914238 
0.00089 
4508.9769 
-89.9 
0.0 



+0.0000003 
0.0000002 
+0.02568 
0.00057 
+0.0001 
0.0001 
+ 1.8 a 
2.7 
b 



0000+ - 0004 
u.uuuu_ 0006 

o ni5+ 0015 

93.2 
110.1 



1 n914239+ - 0000002 
i.uyi4zjy_ .0000004 



n nnns7+ - 02494 

0.0008 /_ .00058 
+0.0002 
0.0001 
+ 180.9 a 
2.7 
b 



4508.9769 
-90.0 
0.0 



-0.0001 



+0.0005 



-0.015 



0.0005 
+0.017 
0.019 

91.4 
108.2 



a These distributions are bimodal with strong peaks 
b By definition.) 

from the Exoplanet Transit database (jPoddanv et al. I 
l2010h and our own secondary eclip se d ata along with 
those of lLopez- Morales et al. I (|2010l) , and lCampo et al.l 
(|2010f ). We exclude the in-transit radial velocity data 
as we do not model fo r the Ross iter-McLaughlin effec t 
(|Gaudi fc Winn ll2007ft. We followlCampo et al. I (l20ll). 
and quote the ILopez- Morales et al. I (|2010f ) eclipse point 
that results from the combined photometry from 1.5 
eclipses, at a single epoch halfway between their observa- 
tions (HJD-2455002.8560 ± 0.0073). T N of course gives 
the transit time to compare to the data, we use e, ujn, 
Tn and P a to calculate the eclipse times, and u>(t) to 
calculate the radial velocity values. We use the MCMC 
techniques explained above to calculate the best-fit pre- 
cessing model, and non-precessing models, except that 
we fit for ecosw and esinw, instead of e and u, as uj is 
poorly constrained as the eccentricity approaches zero. 

We plot our precessing and non-precessing best-fit 
models in Figure [7] and present the MCMC results in 
Table [2] The best-fit models with and without preces- 
sion are similar. The best-fit precessing model features 
a very small rate of precession (Co = 0.003^0.001° d -1 ), 
that barely provides a superior fit once the extra degrees 
of freedom are taken into account (a Bayesian Informa- 
tion Criterion^ of B/C=106.6 for the precessing case, 
compared to £?/C=110.1 for the non-precessing case). 
Thus there is not convincing evidence at this date that 
WASP-12b is precessing. 

Given that th e tim ing offset of the 
ILopez-Morales et al. I (|2010D eclipse detection may 
be suspect, we also refit the non-precessing case with 
this eclipse excluded, and present the MCMC results in 
Table [2] The distribution of eccentricity values fr o m our 
MCMC chain without the ILopez-Morales et al. I (|2010D 
eclipse is non-gaussian (the bottom left panel of Figure[5]) 
and favours a near-zero eccentricity with a tail to higher 
eccentricity values; this limit is e=0.00087io'ooo58- This 

(HJD) of 2455246. 77604±0.00217 (A. Gibson, TRESCA), 
2455253. 32414±0.00287 (F. Lomoz, TRESCA), 2455257.69131 (G. 
Haagen, TRESCA), and 2455265. 33327±0.00129 (H. Kucakova, 

TRESCA). 

13 For the Bayesian Information Criterion (Liddlc 2007) lower- 
values indicate superior fits corrected for the number of free param- 
eters: BIC=x 2 + klnN, where k is the number of free parameters 
and N is the number of data points. 



oj~90° and -90° (where cosoj~0.) 

is due to the fact that although the ecosw values for 
WASP- 12b are well-constrained from the radial- velocity 
data and the combination of the timing of the eclipses 
and transits (the top- left panel of Figure [5]), the esina; 
values are not well-constrained and thus higher eccen- 
tricity values are allowed (the top-right panel of Figure 
[5]) for an argument of periastron where cosw o ^0 at 
oj ~90° and -90° (as can be seen in the contour plot in 
the bottom-right panel of Figure[5]) . Although we are not 
able to rule out higher eccentricity values for WASP- 12b 
with high confidence, the orbit of WASP- 12b is likely 
circular; thus WASP- 12b is no longer an outlier from the 
expectation of the timescale of tidal circularization for 
close-in giant exoplanets. The above analysis would be 
improved by including an a priori constraint on esina; 
using the eclipse durati o n valu es from our own eclipses 
and the iCampo et al. I (120101) Spitzer / IRAC eclipses. 
Unfortunately, although ICampo et al. I ()2010[ ) indicate 
that their best-fit eclipse durations are similar to that of 
the transits and should thus place a t ight constrain on 
esinw near zero, ICampo et al. I ([2010D do not formally 
fit for the duration of the eclipse and do not include the 
associated uncertainties. We discuss the implications of 
fitting our own eclipse durations below. 

4.2. A longer duration secondary eclipse; possible signs 
of material stripped from the planet? 

We also fit our Ks-band photometry and our joint J, 
H & Ks-band photometry with an eclipse model with 
the eclipse duration as a free parameter. Our best-fit 
Ks-band variable eclipse duration fit is presented in Fig- 
ure IH1 Our variable eclipse duration fit does argue for a 
marginally wider secondary eclipse than transit: 

= 1.109+0 039! although this result is only significant at 
the 2.8<7-level. The associated eclipse duration is = 
3.25^0 it hours, longer than the ~2.93 hour optical tran- 
sit found by iHebb et al. I (|2009|) . and l onger than the 
simila r duration IRAC eclipses found by ICampo et al. I 
(2010). That the data suggests a wider secondary eclipse 
than our best-fit model can be seen in the ingress and 
egress of our Ks-band photometry (Figure [3]). 

Our joint analysis of our J, H & Ks-band data also ar- 
gues for a marginally wider secondary eclipse than tran- 
sit: $//// = 1.080;o'°3t or that the duration of the 
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Figure 8. Top-left, top-right and b ottom-left panels: Ma rginalized likelihood for WASP-12b's ecosu, esinu and its eccentricity from the 
non-precessing MCMC chain with the Lopez-Morales ct al. (2010) point excluded. The best-fit value for each panel is given with the solid 
vertical line (for the bottom-left panel this value is nearly indistinguishable from zero), while the 68% credible region is indicated by the 
dotted vertical line. Bottom-right panel: Contour parameter showing the eccentricity, e, and the argument or periastron, uj, of WASP-12b 
again from the same MCMC chain. The 68.3% (1<t; solid-line), 95.5% (2cr; dashed-line), and 99.7% (3cr; dotted-line) credible regions are 
indicated. 

H-band photometry that manifest themselves as a sharp 
decrease in flux, can be well-fit by a significantly wider, 
and deeper secondary eclipse that is unlikely to be phys- 
ical. These facts, combined with a visual inspection of 
Figures I3HS1 suggests that the wider secondary eclipse for 
our joint analysis, is in-fact dominated by our Ks-band 
photometry and the longer duration eclipse may not be 
credible for the joint analysis. 

Our Ks-band photometry is best-fit with a wider sec- 
ondary eclipse than expected: $/j//=1.109toQ39. The 
first possibility to explain this wider than expected 
eclipse is systematic time-correlated, red-noise in our 
photometry, which would not be unexpected as the 
eclipse is only wider than expected at less than the 3<7 
level. Another possibility for this wider eclipse is that the 
planet has a small eccentricity (esinw=0.059lg Q34). We 
have already presented strong evidence that the eccen- 
tric ity of WASP- 1 2 b is q uite likely near zero in H4.ll and 
the lCampo et al.~l (|2010D Spitzer/IRAC eclipse photom- 
etry does not feature a longer duration secondary eclipse. 
Also, although the esinw of the planet is less well con- 
strained in cj4.U an esinw value necessary to explain our 
longer duration eclipse can be ruled out at several sigma 
and thus we find this possibility uncompelling. 

Another possibility - perhaps the most intriguing pos- 
sibility - is that if this apparently wider secondary eclipse 
is not due to systematic effects or due to a small esinw 
for WASP-12b, then it could be due to radiation from 
gas that is escaping from the planet and possibly form- 
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Figure 9. The same as figure [3] except the best-fit model is our 
variable eclipse duration model for our Ks-band photometry. 

eclipse is $// = 3.16^o \q hours. As our J-band data is 
a partial eclipse, it has no ability to constrain the eclipse 
duration on its own. Similarly, as our H-band data suffers 
from significant systematics prior to and during ingress, 
its ability to constrain the eclipse duration is compro- 
mised; in fact, the systematics at the beginning of the 
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ing a cicumstel lar disk. The latter was predicted by 
iLi et al. I (120101). while the fo rmer was arguably recently 
confirmed by iFossati et al. I (|2010ai ) through observa- 
tions that WASP- 12b displayed increased transit depths 
in the UV with COS/HST. An eclipse of this duration 
could argue for material surrounding the planet with a 
projected radius that is approximately 1.9 times the op- 
tical radius of the planet, or at a radius of 3.3 Rj U piter, 
and would thus argue for material emitting radiation 
that is exceeding the Roche lobe, and streaming from 
the planet. This emission could b e due to CO ^2 .292 
bandhead emission, as predicted bv lLi et al.l (|2010D . al- 
though the material around the planet should be cooler 
than the ^4000-5000 K temperatures they predicted for 
the cirumstellar disk and will thus result in reduced emis- 
sion . In the "accret ion stream" hypothesis, advocated 
by (|Lai et al. 1 120101 ) , the material streaming from the 
planet towards the star may be highly localized in a line 
passing through the inner Lagrangian point. The extra 
emission from this stream would be obscured by the star 
earlier than the planet during eclipse ingress and later 
than the planet during eclipse egress. Such a scenario 
is arguably favoured over simply a sphere of evaporating 
material 3.3 Rjupiter in radius, as in the accretion stream 
scenario the emission will arise from a smaller surface 
area; otherwise the Ks-band brightness temperature of 
the planet would have to be anomalously low, given that 
the AFif s =0.309to;°i2% that we observe would have to 
be a combination of the planet and the enveloping mate- 
rial that would have a much larger surface area of emit- 
ting material than the planet itself. 

Alternatively, this wider eclipse could be interpreted as 
the planet passing behind a circumstellar disk that is op- 
tically thick at these wavelengths and extends marginally 
from the star (at least ~1.11 times the stellar radius), 
and therefore obscures the planet earlier and later than 
expected. If the disk is optically thick it will have to 
be due to gas opacity, as the temperature of the disk 
will be well above the dust sublimation temperature, 
and the temperature and density of the disk will have 
to be high enough for the material t o be largely pho- 
toion ized to avoid the "opacity gap" ()Thompson et al. 1 
2005). The disk would also have to be optically thick at 
around 2 /im, but not at the longer wavelengths probed 
by Spitzer/IRAC (3.6 to 8.0 [mi) as the durations of the 
eclips es are not disc r epant from the expected duration 
in the lCampo et al.l (|201Q[ ) photometry. The "accretion 
stream" hypothesis is arguably less contrived, but the 
observed eclipse duration could also result from a com- 
bination of both scenarios. 

An obvious way to differentiate between these two sce- 
narios would be observations of WASP-12 during transit 
in the Ks-band band. If WASP-12 is surrounded by a 
circumstellar disk that is emitting in Ks-band then the 
transit duration will increase. If there is material sur- 
rounding WASP- 12b then its transit will be of the ex- 
pected duration if the material is optically thin, and the 
transit will display an increased depth if the material is 
optically thick. We plan to perform such follow-up ob- 
servations of the transit and eclipse of WASP-12b in the 
Ks and H-bands to differentiate between these various 
scenarios, and to confirm the near-zero eccentricity of 
WASP-12b. Until such follow-up observations take place 




Figure 10. \a (solid-lines), 2a (dashed-lines) , and 3cr (dotted- 
lines) constraints on the Bond albedo and reradiation factor, ftot 
from our Ks, H & J-band secondary eclipse observations of WASP- 
12b. 

we emphasize that our Ks-band photometry is best-fit 
with a wider eclipse at less than the 3<r level. 

4.3. The properties of WASP-12b's atmosphere 

Our measurements of the thermal emission of WASP- 
12b allow us to constrain the characteristics of its atmo- 
sphere, including: its Bond albedo, the level of redistri- 
bution of heat from the day to the nightside at various 
depths, and the planet's dayside bolometric luminosity. 
We parameterize the level of r edistribution by the rera- 
diatio n factor, /, following the lLopez-Morales fc Seagerl 
(2007) definition (i.e. f=j denotes isotropic reradia- 
tion, while f=h denotes redistribution and reradiation 
from the dayside only). Our eclipse depths are con- 
sistent with a range of Bond albedos, A B , and overall 
day to nightside redistribution of heat, f tot (Figure [TU)). 
If we assume a Bond albedo near zero, consistent with 
observations of other hot Jupiters (jCharbonneau et al. I 
Il999t iRowe et alH 120081) and with model predictions 
(|Burrows et al. I I2008D . the best-fit reradiation factor, 
ftot, that results from our three near- infrared eclipse 
measurements is f to t = 0-44lin"o23- This suggests that 
the dayside of WASP-12b reradiates most of the incident 
stellar flux without redistributing it to the nightside. 

As the atmospheres of hot Jupiters may be highly ver- 
tically stratified, different atmospheric layers may redis- 
tribute heat much more or much less efficiently than 
other layers. The best-fit brightness temperatures and 
reradiation factors of the individual atmospheric lay- 
ers probed by our various wavelengths of observations 
are: T B Ks=-2QB&t%K and / A ' s =0.482±°;^° for our Ks- 
band observations, T BH =2765± 7 7 ° 2 K and / K =0.353tg;^ 
for our H-band observations, T B j =2833^173^ an< ^ 
/j=0.389lg Qg 7 for our J-band observations. Our three 
differents bands should be probing high pressure regions, 
deep into the atmosphere of WASP-12b. Specifically if 
the near-infrared opacity is dominated by water vapour 
opacity the J, H fe K-bands should be windows in water 
opacity (jFortnev et al. Il2008l ), and the Ks, H & J-bands 
should be seeing progressively deeper into WASP-12b's 
atmosphere. Within the errors the brightness tempera- 
tures displayed in our three near-infrared bands are simi- 
lar. However, the J and H-band brightness temperatures 
are marginally lower, and taken at face value compared 
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Figure 11. Dayside planet-to-star flux ratios (top) and dayside 
flux at the planet's surface (bottom). The Ks-band (~2.15 fim), ri- 
band (~1.60 fim) and J-band (~1.25 fim) points are our own, whil e 
the z'-band point (~0.9 (an) is from[Lopcz-Moralcs c t al. 1 112010 ). 
Blackbody curves for modest redistribution (/=0.35; T e(J ~2735 K; 
blue dashed line), and dayside only emission (/=!; T eq ~2990 K; 
grey dotted line) are also plotted. We also plot one-dimensio nal, 
radiative transfer spectral models (Fortncy ct al. 2006, 2008) for 
various reradiation factors and with and without TiO/VO. We 
plot models with modest redistribution (/=0.35) with and without 
TiO/VO (magenta-dotted and green-dashed lines, respectively), 
and for dayside only emission (/=!) with and without TiO/VO 
(orange dotted and cyan dot-dashed lines, respectively). The 
models with TiO/VO display temperature inversions. The mod- 
els on the top panel ar e divided by a stellar atmosphere model 
(Hauschildt ct al. 1999) of WASP-12 using the parameters from 
irlebb et al. 112009 ). (Af*=1.35 M e , R*=1.57 R e , T e// =6300 K, 
and log g=4.38). We plot the WIRCam Ks, H and J-band trans- 
mission curves inverted at arbitrary scale at the top of both panels 
(dotted black lines). We integrate our models over the WIRCam 
bandpasses and display the result in the appropriately coloured 
triangles. 



to the Ks-band brightness temperature they suggest a 
modest temperature inversion at very high pressures of 
~100 to 500 mbar, deep in the atmosphere of WASP- 
12b. One explanation for why WASP- 12b might dis- 
play decreased flux at these shorter wavelengths as com- 
pared to the Ks-band, is that the atmospheric depths 
and pressures probed by these shorter wavelength obser- 
vations may be more homogenized than higher altitude 
layers. The efficiency of redistribution of the incident 
stellar flux from the dayside to the nightside should be 
proportional to the ratio of the reradiative (r ro d) to ad- 
vective timescales (T a d v )- It is thought that the reradia- 
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tive timescale should increase with pressure and deptlQ 
The advective timescale is also thought to increase in 
pressure, although it is generally thought that advection 
should win out over reradiation as one des cends through 
the a t mosphere of a typica l hot Jupiter (|Seager et al. I 
120051 iFortnev et al~l l2008f >. Thus, one might expect 
more efficient redistribution of heat at the layers probed 
by our shorter wavelength observations compared to the 
layers probed by our Ks-band observations. Other expla- 
nations for the relatively higher Ks-band emission than 
the H and J-band emission are certainly possible, in- 
cluding: extra flux from a circumstellar disk or mate- 
rial streaming from the planet in the Ks-band, an atmo- 
spheric emission feature at Ks-band, or absorption fea- 
tur es over the H a n d J-b ands. The eclipse depths from 
the iCampo et al. I (|2010l ) Spitzer/IRAC measurements 
will not shed much additional light on this matter, as, if 
water opacity dominates, the Spitzer/IRAC bands to do 
not probe as deeply as the JHK near-infrared bands. 

We compare the depths of our near-infrared eclipses 
to a series of planetary atmosphere models in Figure 111! 
This comparison is made quantitatively as well as qual- 
itatively by integrating the models over the WIRCam 
J, H & Ks band-passes and calculating the x 2 of the 
thermal e mission data compared to th e models. We in- 
clude the lLopez-Morales et al. I (|20100 eclipse depth in 
Figure 111) but do not include it in our \ 2 calculation 
due to the aforemtioned uncertainty with the timing and 
depth of this eclipse . The Spitzer/IRAC eclipse depths 
(|Campo et al. |[2010h are also not included, as of the time 
of writing only the central eclipse times have been re- 
ported. We first plot two blackbody models, the first one 
displaying modestly efficient heat redistribution (/=0.35; 
blue dotted line; T eg ~2735 K), while the latter features 
emission from the dayside only (f=\\ grey dotted line; 

T e(3 ~2990 K). The f=\ blackbody model provides an 
excellent fit to the longer wavelength Ks-band emission, 
and does a reasonable job of fitting our H and J-band 
emission (/=^: X 2 — 5); nevertheless it proves a quanti- 
tatively better fit than the modest redistribution model 
(/=0.35: x 2 =34), which generally underpredicts the ob- 
served emission. 

In Figure [TT] we also compare our measurements to 
a series of one-dimensional, radiative trans fer, spec- 
tral models (|Fortnev et al. I [2001 [2001 [2001 ) with dif- 
ferent reradiation factors that specifically include or 
exclude gaseous TiO/VO into the chemical equilib- 
rium and opacity calculations. In these models when 
TiO/VO are present in gaseous form in the upper at- 
mosphere they act as absorbers at high altitudes and 
lead to hot stratosp heres and temperature inversions 
(|Hubenv et al. 1 120031 ). We present models with mod- 
est redistribution (/=0.35) and dayside only emission 



14 The radiative time-scale (how quickly the planet reradiates 
the incident stellar flux) is thought to be proportional to T ra & ~ 
^gjr ' i* llShowman fc Guillot 1 120021) , where cp is the specific heat 
capacity, a is the Stefan-Boltzmann constant, T is the temperature 
of the atmospheric layer, and g is the gravitational acceleration of 
the planet. 

15 It is thought that the advective timescale (how quickly the 
planet advects the heat to the nightside of the planet; r at j„) can be 
approximated by the radius of the pla net, Rp, divided by the ho r- 
izontal windspeed, U: r adv ~ Rp /U (Showman & Guillot 2002). 



Near-infrared thermal emission from WASP-12b 
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(/=§) wrtn an( i without TiO/VO. The associated x 2 for 
the f=\ models with and without TiO/VO are x 2 =20 
and x 2 =20, while the /=0.35 models with and without 
TiO/VO are x 2 =? 2 and x 2 =83. None of these models 
provide quantitative improvements over the f=\ black- 
body model, as they do not do as good of job of matching 
the longer wavelength Ks-band thermal emission, nor do 
they feature reduced emission in H and J-band. 

Our near-infrared measurements also allow us to es- 
timate the bolometric dayside luminosity of WASP-12b, 
Lday We use a blackbody model with a total reradiation 
factor equal to the best-fit value we calculate from our 
three near-infrared bands (/t O (=0.441); by integrating 
over this model we can estimate Lday as 1.12x10~ 3 Lq. 
Another way of parameterizing the efficiency of the day- 
to-nightside heat redistribution rather than the reradia- 
tion factor is comparing the bolometric dayside luminos- 
ity, Lday, to the nightside luminosity, L nig ht- By follow- 
ing elementary thermal equilibrium calculations one can 
deduce that WASP- 12b should display a total bolometric 
luminosity of L tot = 1.25xlO~ 3 L Q . This suggests that 
89% of the incident stellar irradiation is reradiated by the 
dayside, leaving a mere 11% to be advected to the night- 
side and reradiated. However, caution is encouraged with 
this conclusion as shorter and longer wavelength emission 
for this planet may deviate significantly from that of a 
blackbody. 

4.4. Future Prospects 

We lastly note that the combination of thermal emis- 
sion as prominent as that displayed here with near- 
infrared photometry this precise suggests the possibility 
that thermal phase curve measurements may be possi- 
ble from the ground. For the shortest period exoplanets 
(P~ld or less) even in a single night of observing (8-9 
hours) one could conceivably view the flux maximum of 
the phase curve where hot gas is advected downwind on 
the planet, the decrement in flux during the secondary 
eclipse, and then view a significant fraction of the near- 
sinusoidal decrease as the cool nightside face of the ex- 
oplanet rotates into sight. WASP-12b is an ideal target 
for such observations with its short 1.09d period, and 
its bright dayside emission suggests that thermal phase 
curve observations for this planet should reveal a large 
asymmetry over the course of the orbit as WASP-12b's 
nightside should be cold. Thermal phase curve obser- 
vations from the ground in the near-infrared would re- 
quire one to control the background systematic trends 
that are present in our near-infrared photometry even 
after we correct the flux of our target star with a great 
many reference stars; the feasibility of this task is, as of 
yet, unproven. Nevertheless, we will be investigating the 
possibility of obtaining such near-infrared phase curve in- 
formation in this photometry as well as with future obser- 
vations of WASP-12b. These near-infrared phase curve 
observations will be accompanied by near-simultaneous, 
3.6 and 4.5 /j,m Spitzer/IRAC thermal phase curve ob- 
servations of a full orbit of WASP-12b (P.I. P. Machalek) 
that will allow for an unprecedented understanding of the 
characteristics of the day and nightside deep atmosphere 
of this planet. 

We also plan to reobserve a full, rather than partial, 
eclipse of WASP- 12b in J-band so as to better define 



its thermal emission at that wavelength. Lastly we plan 
to observe the transit of WASP- 12b in the near- infrared 
Ks and H-bands, combined with our aforementioned 
planned reobservations of the eclipse of WASP- 12b in 
these bands. These combined transit and eclipse obser- 
vations will allow us to confirm if the Ks-band eclipse is 
indeed longer in duration than the optical transit, and if 
so whether this is due to material tidally stripped from 
the planet that may or may not form a circumstellar disk 
in this system. 
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